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Abstract. Solar, atmospheric, and reactor neutrino experiments have confirmed 
neutrino oscillations, implying that neutrinos have non-zero mass, but without pinning 
down their absolute masses. While it is established that the effect of neutrinos on 
the evolution of cosmic structure is small, the upper limits derived from large-scale 
structure could help significantly to constrain the absolute scale of the neutrino masses. 
In a recent paper the 2dF Galaxy Redshift Survey (2dFGRS) team provided an upper 
limit m^tot < 2.2 eV , i.e. approximately 0.7 eV for each of the three neutrino flavours, 
or phrased in terms of their contribution to the matter density, f2„/f2 m < 0.16. Here 
we discuss this analysis in greater detail, considering issues of assumed 'priors' like the 
matter density f2 m and the bias of the galaxy distribution with respect to the dark 
matter distribution. As the suppression of the power spectrum depends on the ratio 
f2„/f2 m , we find that the out-of-fashion Mixed Dark Matter model, with f2„ = 0.2, 
fl m = 1 and no cosmological constant, fits both the 2dFGRS power spectrum and the 
CMB data reasonably well, but only for a Hubble constant Hq < 50 kms _1 Mpc _1 . 
As a consequence, excluding low values of the Hubble constant, e.g. with the HST 
Key Project, is important in order to get a strong upper limit on the neutrino masses. 
We also comment on the improved limit obtained by the WMAP team, and point out 
that the main neutrino signature comes from the 2dFGRS and the Lyman a forest. 
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1. Introduction 

The wealth of new data from e.g. the cosmic microwave background (CMB) and large- 
scale structure (LSS) in the last few years indicate that we live in a flat Universe 
where ~ 70 % of the mass-energy density is in the form of dark energy, with matter 
making up the remaining 30 % . The WMAP data combined with other large-scale 
structure data [Tj gives impressive support to this picture. Furthermore, the baryons 
contribute only a fraction /b = f2b/fi m ~ 0.15 (fib and f2 m are, respectively, the 
contribution of baryons and of all matter to the total density in units of the critical 
density p c = 3H%/8nG = 1.879 x 10~ 29 /i 2 gem" 3 , where H = lOO/i kms^Mpc^ 1 is 
the present value of the Hubble parameter) of this, so that most of the matter is dark. 
The exact nature of the dark matter in the Universe is still unknown. Relic neutrinos 
are abundant in the Universe, and from the observations of oscillations of solar and 
atmospheric neutrinos we know that neutrinos have a mass El IH 13 13 E] and 
will make up a fraction of the dark matter. However, the oscillation experiments can 
only measure differences in the squared masses of the neutrinos, and not the absolute 
mass scale, so they cannot tell us how much of the dark matter is in neutrinos. From 
general arguments on structure formation in the Universe we know that most of the 
dark matter has to be cold, i.e. non-relativistic when it decoupled from the thermal 
background. Neutrinos with masses on the eV scale or below will be a hot component 
of the dark matter. If they were the dominant dark- matter component, structure in 
the Universe would have formed first at large scales, and smaller structures would form 
by fragmentation (the 'top-down' scenario). However, the combined observational and 
theoretical knowledge about large-scale structure gives strong evidence for the 'bottom- 
up' picture of structure formation, i.e. structure formed first at small scales. Hence, 
neutrinos cannot make up all of the dark matter (see for a review). Neutrino 
experiments give some constraints on how much of the dark matter can be in the 
form of neutrinos. Studies of the energy spectrum in tritium decay ^U] provide an 
upper limit on the electron neutrino mass of 2.2 eV (95 % confidence limit). For the 
effective neutrino mass scale involved in neutrinoless double beta decay an upper limit 
of 0.34 eV (90 % confidence) has been inferred [TT] but then under the assumptions that 
neutrinos are Majorana particles (i.e. their own antiparticles), and the translation from 
the effective neutrino mass scale in neutrinoless double beta decay to neutrino mass 
eigenvalues requires assumptions about the neutrino mass hierarchy and the CP phases 
in the neutrino mixing matrix. 

From cosmology, an analysis the 2dF Galaxy Redshift Survey [12] found m u ^ot < 
2.2 eV as an upper limit on the sum of the (degenerate) mass eigenvalues. In 
[TB*] the 2dFGRS was combined with pre- WMAP CMB data to give an upper limit 
m^tot < 0.9 eV. The WMAP team P improved this result to m u ^ tot < 0.71 eV (95 
% confidence) from a combination of WMAP, ACBAR [H], CBI [TS] for the CMB, 
and 2dFGRS and the power spectrum inferred from Lyman a forest [To^ ITTj for large 
scale structure, a factor of roughly three better than the 2dFGRS limit (not an order 
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of magnitude as stated in [I]). This limit is comparable to what the pioneering study 
in [T5] predicted would be possible with the Sloan Digital Sky Survey [TH] and WMAP 
[TH] . and have implications for neutrino oscillation experiments as it seems to call into 
question the Liquid Scintillator Neutrino Detector (LSND) result [20], where the mass- 
square difference involved was ~ 1 eV 2 [2HE2]- (However, as pointed out in [2H] it 
is premature to say that cosmology rules out the LSND results.) Note that neutrinos 
with eV masses are basically indistinguishable from cold dark matter at the epoch of 
last scattering, and therefore they have little effect on the CMB. The important role of 
the WMAP data in the cosmological neutrino mass limit is to break degeneracies in the 
parameter space that will otherwise limit the ability to constrain neutrino masses from 
the large-scale structure data. 

In this paper we discuss in detail the recent cosmological neutrino mass limits, 
concentrating on the 2dFGRS and the WMAP + 2dFGRS limits and the various 
parameter degeneracies involved in the analysis. In particular we discuss the role of the 
bias of the galaxy distribution with respect to the mass distribution, non-linear effects, 
and the necessity of using independent information about cosmological parameters like 
Q m and h ('priors'). We will throughout this paper work within the context of flat ACDM 
models, which are favoured by a wealth of observational data PjJ [2H Q3|, however we 
also comment on Mixed Dark Matter (MDM) models in discussion of the analysis of the 
2dFGRS data. In fact we find that an MDM model can still provide a reasonable fit to 
the 2dFGRS and WMAP data, although with a low Hubble constant, so that external 
constraints on the Hubble constant are important in order to get a strong upper limit 
on the neutrino masses. 

The structure of this paper is as follows: in section 2 we give a brief overview of how 
neutrinos affect structure formation in the Universe. In section 3 we consider galaxy 
redshift surveys as a probe of neutrino masses, starting with a brief summary of the 
analysis in [12 . Since the exact relationship between the distribution of the galaxies and 
that of the dark matter is unknown, we discuss different ways of taking this uncertainty 
into account. We also discuss the role of priors on parameters degenerate with massive 
neutrinos. In section 4 we give a brief overview of other cosmological probes of neutrino 
masses before we summarize and conclude in section 5. 



2. Massive neutrinos and structure formation 



The relic abundance of neutrinos in the Universe today is straightforwardly found 
from the fact that they continue to follow the Fermi-Dirac distribution after freeze- 
out, and their temperature is given in terms of the CMB temperature Tcmb today as 

% = (VII^Tcmb, 



n - ^T 3 (1) 

n v — 2 CMB) \ l ) 



where £(3) ~ 1.202, which gives n v w 112 cm 3 at present. Neutrinos are so light that 
they were ultra-relativistic at freeze-out. Their present contribution to the mass density 
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can therefore be found by multiplying n v with the total mass of the neutrinos m„ )tot , 
giving 

Qh 2 - mu ' tot (2) 
v ~ 94 eV 1 j 

for Tcmb = 2.728 K. Several effects could modify this simple relation. If any of the 

neutrino chemical potentials were initially non-zero, or there were a sizable neutrino- 

antineutrino asymmetry, this would increase the energy density in neutrinos and give 

an additional contribution to the radiation energy density. However, from Big Bang 

Nucleosynthesis (BBN) one gets a very tight limit on the electron neutrino chemical 

potential, since the electron neutrino is directly involved in the processes that set the 

neutron-to-proton ratio. Also, within the standard three-neutrino framework one can 

extend this limit to the other flavours as well. The recent results of the KamLAND 

experiment [2EJ confirmed the Large Mixing Angle (LM A) solution for the solar neutrino 

oscillations, and combined with the atmospheric data indicating maximal mixing in 

this sector, it has been shown that flavour equilibrium is established between all three 

neutrino species before the epoch of BBN |26ll2T| l2*%]. so that the BBN constraint on the 

electron neutrino asymmetry applies to all flavours, which in turn implies that the lepton 

asymmetry cannot be large enough to give a significant contribution to the radiation 

energy density. Recent analyses of WMAP and 2dFGRS data give independent, 

although not quite as strong, evidence for small lepton asymmetries |2*3|I29|. Within the 

standard picture, equation (JTJ) should be accurate, and therefore any constraint on the 

cosmic mass density of neutrinos should translate straightforwardly into a constraint on 

the total neutrino mass, according to equation (J2J). If a fourth, light 'sterile' neutrino 

exists, sterile-active oscillations would modify this conclusion. No sterile neutrinos are 

required to explain the solar and atmospheric neutrino oscillation data jHO] , and the only 

hint so far comes from the possible detection of — » V e oscillations with a small mixing 

angle and a mass- square difference ~ 1 eV 2 at the LSND |Uj. Since there are only two 

independent mass-squared differences in the standard three-neutrino scenario, and they 

are orders of magnitude smaller, this hints at the existence of a fourth, light sterile 

neutrino. However, as said, this has little support in the solar and atmospheric data. 

The status of the LSND results will in the near future be clarified by the MiniBooNE 

experiment |3*T] . 

Finally, we assume that the neutrinos are nearly degenerate in mass. Current 
cosmological observations are sensitive to neutrino masses ~ 1 eV or greater. Since 
the mass-square differences are small, the assumption of a degenerate mass hierarchy 
is therefore justified. This is illustrated in figure ^ where we have plotted the mass 
eigenvalues mi,m2, as functions of m U) tot = mi + m2 + m% for Am^ = 5 x 10~ 5 eV 2 
(solar) and Am^ = 3 x 10~ 3 eV 2 (atmospheric), for the cases of a normal hierarchy 
(mi < ui2 < 1TL3), and an inverted hierarchy (m^ < m\ < 1712). As seen in the figure, for 
m u,tot > 0.4 eV the mass eigenvalues are essentially degenerate. 

We will in this paper look at cosmological models with four components: baryons, 
cold dark matter, massive neutrinos, and a cosmological constant. Furthermore, we 
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Figure 1. Neutrino mass eigenvalues as functions of their total for the cases 
of normal (top panel) and inverted (bottom panel) hierarchies. The vertical line 
marked 'oscillations' is the lower limit derived from the measured mass-squared 
differences for the two hierarchies. The other vertical lines are upper limits from 
WMAP+CBI+ACBAR+2dFGRS+Ly a, 2dFGRS, and 3 H j3 decay. 

restrict ourselves to adiabatic, linear perturbations. The basic physics is then fairly 
simple. A perturbation mode of a given wavelength A can grow if it is greater than 
the Jeans wavelength Aj determined by the balance of gravitation and pressure, or rms 
velocity in the case of massless particles. Above the Jeans scale, perturbations grow 
at the same rate independently of the scale. Long after matter-radiation equality, all 
interesting scales are above Aj and grow at the same rate, and in models where all 
the dark matter is cold, the time and scale dependence of the power spectrum can 
therefore be separated at low redshifts. Light, massive neutrinos can, however, move 
unhindered out regions below a certain limiting length scale, and will therefore tend 
to damp a density perturbation at a rate which depends on their rms velocity. The 
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presence of massive neutrinos therefore introduces a new length scale, given by the size 
of the co-moving Jeans length when the neutrinos became non-relativistic. In terms of 
the comoving wavenumber, this is given by 



for three equal-mass neutrinos. The growth of Fourier modes with k > k ni will be 
suppressed because of neutrino free-streaming. The free-streaming scale varies with the 
cosmological epoch, and the scale and time dependence of the power spectrum cannot 
be separated, in contrast to the situation for models with cold dark matter only. 

The transfer functions of the perturbations in the various components provide a 
convenient way of describing their evolution on different scales. Using the redshift z to 
measure time, the transfer function is formally defined as 



where S(k, z) is the density perturbation with wavenumber k at redshift z, and D is the 
linear growth factor. The normalization redshift z* corresponds to a time long before 
the scales of interested have entered the horizon. The transfer function thus gives the 
amplitude of a given mode k at redshift z relative to its initial value, and is normalized 
so that T{k — 0, z) — 1. The power spectrum of the matter fluctuations can be written 
as 



where P*(k) is the primordial spectrum of matter fluctuations, commonly assumed to 
be a simple power law P*(k) = Ak n , where A is the amplitude and the spectral index 
n is close to 1. It is also common to define power spectra for each component, see [22] 
for a discussion. Note that the transfer functions and power spectra are independent of 
the value of the cosmological constant as long as it does not shift the epoch of matter- 
radiation equality significantly. 

Accurate determination of the transfer function requires the solution of the coupled 
fluid and Boltzmann equations for the various components. This can be done using 
one of the publicly available codes, e.g. CMBFAST [33] or CAMB [3J]. Analytical 
approximations are also available, and they are very useful when one wants very quick 
computation of transfer functions. Accurate fitting formulas for the transfer function 
were derived by [32]. These analytic approximations are good at realistic baryon 
fractions, i.e. 0.1-0.2, with the errors typically smaller than 4 %. In figure El we show 
the transfer functions for models with Q m = 0.3, J\ = 0.04, h = 0.7 held constant, but 
with varying neutrino fraction. One can clearly see that the small-scale suppression of 
power becomes more pronounced as the neutrino fraction f v = VLy/Vl^ increases. The 
suppression of the power spectrum on small scales is roughly proportional to f u : 




(3) 




(4) 



P m (k,z) = P*(k)T 2 (k,z), 



(5) 



AP m (AQ 

Pra(k) 



8fu- 



(6) 
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Figure 2. Ratio of the transfer functions for various values of f2„ to the one for 
£l v = 0. The other parameters are fixed at O m = 0.3, f?b = 0.04, h = 0.7. 

3. Constraining the total neutrino mass with the 2dFGRS 

In an earlier short paper ^2] we used the power spectrum of the galaxies as measured 
by the 2dFGRS to limit the fractional contribution f v = Q u /Q m to the matter density 
of massive neutrinos, and on their total mass m Uttot = 9AQ u h 2 eV. The present section 
starts with a short summary of the 2dFGRS and the analysis in jT2] before going into 
a more detailed discussion of the various ingredients involved in the analysis. In pQ 
the WMAP team derived a stronger limit on the total neutrino mass than what was 
obtained from the 2dFGRS + various priors. However, the 2dFGRS power spectrum 
played a central role in the WMAP neutrino mass limit. As the CMB is insensitive 
to neutrino masses in the eV range, the main role of the WMAP data is to provide 
tight constraints on parameters that may otherwise partly mimic the effect of massive 
neutrinos on the matter power spectrum. Therefore our discussion of priors should also 
be of interest in understanding how the WMAP limit was obtained. 

3.1. The 2dF Galaxy Redshift Survey 

The 2dF Galaxy Redshift Survey [35^ has measured the redshifts of more than 230 000 
galaxies with a median redshift of z m w 0.11. One of the main goals of the survey was 
to measure the galaxy power spectrum on scales up to a few hundred Mpc, thus filling 
in the gap between the small scales covered by earlier galaxy surveys and the largest 
scales where the power spectrum is constrained by observations of the CMB. A sample 
of the size of the 2dFGRS survey allows large-scale structure statistics to be measured 
with very small random errors. An initial estimate of the convolved, redshift-space 
power spectrum of the 2dFGRS has been determined jSH] based on a sample of 140 
000 redshifts. On scales 0.02 < k < 0.15h Mpc -1 the data are robust and the shape 
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Figure 3. Ratio of power spectra for Vt y = 0.01 (bottom line) and Q, u — 0.05 (top line) 
to the one for f2„ = (horizontal line) with amplitudes fitted to the 2dFGRS power 
spectrum data (vertical bars) in redshift space. We have fixed f2 m = 0.3, h = 0.7, and 
ilbh 2 = 0.02. The vertical dashed lines limit the range in A; used in the fits. 



of the power spectrum is not affected by redshift-space or nonlinear effects, though the 
amplitude is increased by redshift-space distortions. One should bear in mind that 
what is measured is the convolution of the true galaxy power spectrum with the window 
function of the survey [31H I3T|. 



where W is the window function. 

As an illustration of the potential of the 2dFGRS to constrain neutrino masses, 
we show in figure El the ratio of the power spectra for Q u = 0.01, and 0.05 (all other 
parameters are fixed at their 'concordance model' values given in the figure caption) to 
the power spectrum for Q u = after they have been convolved with the survey window 
function, and their amplitudes fitted to the 2dFGRS power spectrum data. It is seen 
from figure H3 that the error bars on the power spectrum data points are correlated, 
as discussed in [36J, and this is taken into account throughout this paper by using the 
full covariance matrix of the data when computing likelihoods. For the 32 data points, 
the VL V = 0-model had x 2 — 32.9, Q u = 0.01 gives \ 2 — 33.4,whereas the model with 
Vt v = 0.05 provides a poor fit to the data with \ 2 — 92.2. 

3.2. Previous results from the 2dFGRS 

In ^2] six parameters were used to describe the matter power spectrum: 

• The fraction of massive neutrinos (hot dark matter) f v = Q„/Q m . 

• The combination Q m h, which describes the shape of the cold dark matter power 
spectrum. 




(7) 
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• The baryon fraction f^ = fl^/Q^. 

• The present value of the Hubble parameter h. 

• The scalar spectral index of the primordial density perturbation spectrum, n. 

• The amplitude A of the galaxy power spectrum. 

The amplitude is a free parameter to take into account the fact that what is measured is 
the power spectrum of galaxies, not of all the matter, and that it is measured in redshift 
space, not in real space. The latter effect has been shown in to correspond to a shift 
in the amplitude of the spectrum. The first effect is parametrized by the so-called bias 
parameter 

The scale-dependence of the bias factor is not well known. We will discuss this in 
more detail in a later subsection, for the time being we note that one would expect 
the relation between the distribution of dark matter and luminous matter to be simple 
on large scales. This is borne out by numerical simulations [38J, and two independent 
analyses have shown that the 2dFGRS power spectrum is consistent with a constant 
bias on the scales relevant for our analysis HH] . Thus, we took the redshift space 
distortions and the bias into account by leaving the amplitude A of the power spectrum 
as a free parameter, which means that we used the shape of the 2dFGRS power spectrum, 
and not its amplitude, to constrain f u . 

The main effect of the massive neutrinos is to reduce the power on scales smaller 
than the neutrino free-streaming scale. This effect may, however, be partially masked 
by other effects. Obviously, lowering the amplitude or the scalar spectral index n will 
reduce the power. Also, the baryon fraction and Q m h interfere with the neutrino signal. 
Therefore, constraints on the neutrino mass from the galaxy power spectrum depends on 
the information we have about other parameters ('priors'). In ^2] we added constraints 
from independent cosmological probes: a Gaussian prior on the Hubble parameter 
h = 0.70 ± 0.07, consistent with the results from the HST Hubble Key Project |41j . 
a Gaussian prior Q^h 2 = 0.020 ± 0.002 on the physical baryon density from Big Bang 
Nucleosynthesis (321) an d a Gaussian prior n = 1.0 ±0.1 on the scalar spectral index. 
Furthermore, we considered two different priors on fl m : 

• Under the assumption of fl m + = 1, a Gaussian prior Q m = 0.28 ± 0.14 was 
obtained from surveys of high redshift Type la supernovae |4*31 14*4"] . 

• A uniform ('top hat') prior in the range 0.1 < Q m < 0.5. Given our prior on h, 
Q m < 0.5 ensured that the ages of the Universes in the models considered were 
greater than 12 Gyr. 

As noted earlier, the transfer function does not depend on Q\ and so the assumption 
Q m + Qj^ = 1 enters only through the Supernova Type la prior on Q m . 

For each set of parameters, we computed the theoretical matter power spectrum, 
and obtained the y 2 for the model given the 2dFGRS power spectrum. We then 
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Figure 4. Probability distributions, normalized so that the area under each curve is 
equal to one, for f v with marginalization over the other parameters, as explained in 
the text, for N v = 3 massive neutrinos and f2 m /i = 0.16 (full line), 0.20 (dotted line), 
and 0.24 (dashed line). 

calculated the joint probability distribution function for f v and V = Q m h (which 
represents the shape of the CDM power spectrum) by marginalizing over A, h and /b 
weighted by the priors given above. For A we used a uniform prior in the interval 
0.5 < A < 10, where A = 1 corresponds to the normalization of the 'concordance 
model', discussed in jHOj- Using instead a prior uniform in log A, or fixing A at the 
best-fit value had virtually no effect on the results. Figure 0] shows the probability 
distributions for f v for three different values of Q m h. Marginalizing over Q m h using the 
uniform prior on Q m , we got an upper limit f v < 0.13 at 95% confidence for n = 1. 
Increasing n increases power on small scales and leaves more room for suppression by 
the massive neutrinos, and and upon marginalizing over the full range of n with a prior 
n = 1.0 ± 0.1 we found f v < 0.16 at 95 % confidence. For Q m h 2 = 0.15 this corresponds 
to a total neutrino mass m Uttot = 2.2 eV. The results with the Supernova Type la 
prior on fl m were identical. For comparison, marginalizing without any priors, the limit 
becomes f v < 0.24. Adding just a prior on Q m , we find f v < 0.15, so this is clearly 
the most important prior. Marginalizing with just a prior on h or on Q^h 2 , the 95 % 
confidence limit becomes f v < 0.20. Clearly the priors play a crucial role in the analysis, 
and we will discuss their role in the next subsections. To facilitate the comparison with 
the WMAP analysis in pQ we will from now on carry out our analysis in terms of the 
physical densities Ui = Qih 2 , where i = m, u, b. 

3. 3. The prior on u m 

As noted above, the prior on the matter density is crucial, and the tight correlation 
between m Uttot and u m is illustrated in figure El In fact, without a prior on u m no non- 
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Figure 5. Confidence contours (68 and 95 %) from the 2dFGRS data alone in the 
plane of total neutrino mass TO„ jt ot = 94u>„ eV and the physical matter density u> m . 
The bias parameter and as have been marginalized over with top-hat priors, tub, h, 
and n are fixed at their best-fitting values. 

trivial upper limit on uj v is obtained (but note that one still finds an upper limit on f u ). 
WMAP provides a constraint uo m = 0.14 ±0.02 for spatially flat ACDM models, but it is 
interesting to note that one from the CMB and 2dFGRS alone cannot rule out models 
with Q m = ficDM + fib + fi^ = 1. To illustrate this point, we consider the following 
three models, all with uuh = 0.024: 

(i) A Mixed Dark Matter (MDM) model with fi m = 1, Vt u = 0.2, h = 0.45, and 
n = 0.95. The neutrino mass fraction is thus /„ = 0.2. 

(ii) A ACDM model with fi m = 0.3, Vt u = 0, h = 0.7, and n = 1.0. 

(iii) A pure CDM model with tt m = l,Q u = 0,h = 0.45, and n = 0.95. 

For the pre- WMAP CMB, we use the recent compilation in [13] • Model 1 has 
X 2 /datum =1.14 for the CMB (28 points) and X 2 /datum = 1.13 for the 32 2dFGRS 
power spectrum data points. Model 2 has x 2 /datum = 1.11 for the CMB, x 2 — 1-03 for 
the 2dFGRS. The WMAP data discriminate better between these two models, Model 
1 having % 2 /datum = 1.15 (899 points) and model 2 % 2 /datum = 1.08, but we see 
from figure El that the models are look reasonable and note that we have not carried 
out any systematic search for a best-fitting model of the three types. Thus, the first 
two models seem to offer acceptable descriptions of the CMB and 2dFGRS data, and 
from these data alone MDM is still a viable alternative to the 'concordance' ACDM. 
So why did several pre- WMAP studies find that the CMB and 2dFGRS prefer a low 
matter density and a cosmological constant ? This is because they considered (very 
reasonably) neutrinos to be essentially massless. Model 3 illustrates the point: it is a 
pure CDM model with massless neutrinos. It gives a reasonable description of the CMB 
data, but has \ 2 — 67 for the 2dFGRS data points, and hence it will be disfavoured in a 
joint analysis. The CMB cannot distinguish between eV-mass neutrinos and cold dark 
matter, and hence model 1 and 3 provide comparable descriptions of the CMB data, 
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but the galaxy power spectrum does distinguish between the two, and massive neutrinos 
provides the necessary reduction in small-scale power to fit the data points. Of course, 
one cannot look at the CMB and 2dFGRS data alone, and it is not our intention to 
'resuurect' the MDM model: it has problems with the evolution of cluster abundances 
with redshift [IE], needs a low value of the Hubble parameter, and Q m = 1 is clearly at 
variance with independent measurements from e.g. the baryon fraction in clusters |47j . 
but we wish to make the two following points: 

• The CMB + 2dFGRS alone cannot rule out MDM; one needs a further prior, e.g. 
the Hubble constant from HST j¥T], supernovae Type la [HU |Hj, a prior on f2 m 
from the cluster baryon fraction or evolution of cluster abundance with redshift 

m- 

• The statement that CMB + 2dFGRS alone provides evidence independent of that 
of supernovae Type la for a cosmological constant / dark energy is too strong. This 
result is obtained only when neutrino masses are assumed to be negligible. If one 
allows for massive neutrinos, acceptable fits to the CMB+2dFGRS can be obtained 
with Q m = 1 and Q u ~ 0.2. 

3.4- The prior on the Hubble parameter 

We saw in the previous subsection that if we low values of the Hubble constant, h < 0.5, 
MDM models provide reasonable fits to the CMB and 2dFGRS power spectra, but by 
combining with the HST prior on h, one will obtain the by now usual 'concordance' 
values Q m ~ 0.3, ~ 0.7. The derivation of a strong upper limit on the total neutrino 
mass therefore depends on our ability to exclude values of h much below 0.7. Figure [3 
shows the degeneracy between h and u u . If we allow f2 m = 1, and drop the prior on h, 
we find that a universe with a low value of h < 0.5 is a viable option if the total neutrino 
mass is a few eVs. Note that these models also have ages > 12 Gyr, consistent with 
ages of globular clusters jlHj. Thus, without a prior on h, the CMB prior on f2 m would 
have been weaker, and which would also have affected the upper limit on the neutrino 
masses. 

3. 5. The prior on the scalar spectral index 

As noted earlier, there is also a degeneracy between the scalar spectral index n and 
u u , illustrated in figure |H1 Motivated by the pre- WMAP CMB data we used a prior 
n = 1.0 ± 0.1 in [TJ. For flat ACDM models the WMAP data give n = 0.99 ± 0.04 
(see table 1 in [T]), but in their full analysis including other datasets there is some 
evidence for a running scalar spectral index. However, it has been argued that this may 
be because of their treatment of the Lyman a forest power spectrum, and that in a 
more conservative approach one finds that a scale-invariant primordial power spectrum, 
n = 1 is consistent with the data jUJj. We will therefore not consider a running spectral 
index in this paper. 
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Figure 6. MDM (n m = 1, Cl v = 0.2, h = 0.45, n = 0.95), CDM (Q m = 1, Q, v = 0, 
h = 0.45, ri = 0.95) and ACDM (Q m = 1, fi A = 0.7, fl v = 0, h = 0.7, n = 1.0) models 
(all with f2b/i 2 = 0.024) compared with data. The upper panel shows the pre- WMAP 
data, the middle panel the WMAP data and the lower panel the 2dFGRS data along 
with the three models considered. We have normalized the models to each data set 
separately, but otherwise these are assumed models, not formal best fits. 
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Figure 7. Confidence contours (68 and 95 %) from the 2dFGRS data alone in the 
rnv,tot~ h plane. The bias parameter and as have been marginalized over with top-hat 
priors, while u) m , wt,, and n have been fixed at their best-fitting values. 




Figure 8. Confidence contours (68 and 95 %) from the 2dFGRS data alone in the 
TO v,tot -n plane. The bias parameter and erg have been marginalized over with top-hat 
priors, while u) m , Uh and h have been fixed at their best-fitting values. 

3. 6. The prior on the baryon density 

There is also some degeneracy between the baryon density and the neutrino masses, 
as shown in figure 01 because increasing the baryon density suppresses power on small 
scales. The degeneracy is, however, less serious than for the other parameters, and 
WMAP provides a tight constraint on from the ratio of the amplitudes of the first 
and second peaks in the CMB power spectrum, = 0.024 ± 0.001, consistent with 
standard BBN g2]. 
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Figure 9. Confidence contours (68 and 95 %) from the 2dFGRS data alone in the 
m^tot - Wb plane. The bias parameter and as have been marginalized over with top-hat 
priors, while w m , h, and n have been fixed at their best-fitting values. 



3. 7. Non-linear fluctuations 

At small scales, one eventually enters the regime where the power spectrum is no longer 
linear. Therefore, as a further test of the stability of our analysis in ^2] , we used the 
full set of priors, but only the power spectrum data at scales k < 0.1 /iMpc -1 and 
found that the limit increased to /„ < 0.20. Cutting the power spectrum at this scale 
is, however, very conservative. To see this, we follow the analysis of non-linearities in 
[HH| . Defining 

A 2 (^) = ^P(fc), (9) 

the crossover from linear to non-linear behaviour is at the co-moving momentum k cut 
where A 2 (fc cut ) = 1. This corresponds roughly to the point where A^ on _ linear — Af inear ~ 
A 2 inear . We use the approximate relation between the the linear and non-linear spectrum 
found by [ST]. The cut in the spectrum should be made where the non-linear effects 
are of the same order as the suppression of small scale-power from massive neutrinos, 
which is given by equation (JHJ). For neutrino masses ~ 1 eV, this gives A 2 inear ~ 3. We 
chose k cut = 0.15 /iMpc -1 in our analysis, and at that point A 2 inear « 0.7. Crude as this 
argument may be, it clearly indicates that the neutrino mass signal dominates possible 
non-linear effects in the 2dFGRS power spectrum data used in our analysis. We note that 
the linear matter power spectrum is convolved with the survey window function before 
comparison with the data. However, for values of k greater than ~ 0.1 /iMpc -1 the 
window function is sharply peaked at k, and so there is little mixing with smaller scales 
[3*7] . The WMAP team went one step further and considered non-linear corrections at 
the smallest scales in the analysis. Nevertheless, one would expect these effects to be 
small for k < 0.15 /iMpc -1 . 
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3.8. Scale- dependent bias 

The power spectrum of the galaxy distribution might be 'biased' relative to the matter 
power spectrum, hence it might introduce systematic error in the estimation of the 
neutrino mass. 

Indeed it is well established that on scales less than ~ 10 Mpc different galaxy 
populations exhibit different clustering amplitudes, the so-called morphology-density 
relation (e.g. |H21 E31 E3J ESj)- Hierarchical merging scenarios also suggest a more 
complicated picture of biasing as it could be non-linear, scale-dependent and stochastic 
(e.g. |2HlEniEZlEHlEniinniini]- But is the biasing still scale-dependent at the large scales 
{k~ l > 7 h~ l Mpc) where we analyse the 2dFGRS power spectrum ? Let us consider 
the ratio of galaxy to matter power spectra, and use the ratio of these to define the bias 
parameter as in equation (jHJ). To illustrate the dramatic effect that scale-dependent 
might have we assume the following simple form: 

b(k) = alog(k/k*) + c, (10) 

where we fix k* = 0.15 /iMpc -1 (note that a shift in k* can be absorbed in a change in 
c), but allow a and c to vary. 

Analysis of the semi-analytic galaxy formation models in shows that on large 
scales the biasing function b(k) is nearly constant to high degree. In our parameterization 
()10|) even the brightest galaxies (L > 0.75L*, where is the characteristic luminosity 
of the Schechter luminosity function) are weakly biased, with slope a < 0.15 over the 
scales 3 < k^ 1 < 16/i _1 Mpc [02] • The simulations in ;60] also suggest scale-independent 
biasing on scales larger than 10 h -1 Mpc at late times. 

Observationally, the bi-spectrum analysis of the 2dFGRS showed almost no 
deviation from linear biasing jlOj and combined analysis of 2dFGRS with CMB data on 
scales of 0.02 < k < 0.15 h Mpc -1 jSHJ gave b ~ 1 for L* galaxies. Furthermore, the ratio 
of the power-spectra of blue and red galaxies in 2dFGRS [B^ is almost constant over the 
range of our analysis, 0.02 < k < 0.15 /iMpc -1 . This suggests (as a necessary, but not 
sufficient condition) that the galaxy power spectrum is proportional to the underlying 
matter power spectrum. Based on these theoretical and observational studies we argue 
that scale-dependent biasing is unlikely to pose a problem in estimating the neutrino 
mass from the 2dFGRS. 

We then redo the analysis of the 2dFGRS P g {k) with u u , a and c as free parameters. 
The remaining parameters are fixed at their 'concordance' values (Q m = 0.3, h = 0.7 
etc.). We distinguish between two cases: a < (bias increasing with length-scale and 
a > (bias decreasing with length-scale). In the first case, the best fit is found at 
uj v = 0, a = 0, whereas in the second case one finds the best fit for uj v « 0.029, a = 0.72, 
c = 1.1, i.e. for a non-zero neutrino mass m^ tot ~ 2.8 eV. This is understandable, since 
b{k) in this case has the opposite effect of massive neutrinos, so the two effects can be 
'tuned' to give a very good fit to the 2dFGRS data with a non-zero neutrino mass which 
is unrealistically high. Figure ITU1 shows the likelihood contours in the m u tot -a plane. As 
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Figure 10. Likelihood contours (68, 95, and 99 %) in the m Uj t t-a plane after 
marginalizing over c, where a and c are defined in equation (|10fl . The horizontal 
line a = 0.15 is an upper limit estimated from simulations of biasing (see text). 

argued above, simulations and observations argue for a < 0.15 at large scales, and in 
what follows we shall assume constant biasing. 

3.9. The case of scale-independent bias in more detail 

WMAP provides tight constraints on u m , h, n, and u\,, and we have seen in previous 
subsections that having good constraints on these parameters is essential to obtaining a 
good upper limit on m^tot- Also, from the CMB one can constrain the amplitude of the 
matter power spectrum (quoted, e.g. in terms of the rms mass fluctuations in spheres 
of radius 8 h^ 1 Mpc, erg). If the biasing and redshift-space distortions were known, then 
this would translate directly to a constraint on the amplitude of the 2dFGRS power 
spectrum, and could potentially tighten the constraint on m u tot . The WMAP analysis 
makes use of the constraint on the amplitude of the matter power spectrum from the 
CMB by introducing a prior on the bias parameter, taking it to be a Gaussian with 
b = 1.04 ±0.11 [HE] as found in the analysis of the 2dFGRS bispectrum [10] • As the 
analysis in j^U] was performed for a different range of scales than those involved in the 
analysis of the linear part of the 2dFGRS power spectrum, and did not take neutrino 
masses into account (even though the cosmology dependence in the bispectrum analysis 
is mild) and questions have been raised about this approach [HS], it is worthwhile to 
take a closer look at the WMAP approach. To do this, we need to go into the issue of a 
constant bias and redshift-space distortions in more detail than in the previous sections, 
and we will do so following [35] . The WMAP analysis was more detailed [HH], but we 
think that our simplified version captures the main points. 

We now need to take two effects into account: the fact that the 2dFGRS power 
spectrum is given in redshift space, and that the galaxy distribution may be biased with 
respect to the mass distribution. At redshift z = 0, the relation between the real-space 
normalization at redshift 0, erg m (0), of the matter power spectrum and the normalization 
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of the galaxy power spectrum is given by 

a 8 R g (L s ,0) = &(L a ,0)<r 8m (0), (11) 

where L s « 1.9L* for the 2dFGRS. We follow jSlj and assume that galaxy clustering 
evolves weakly in the range of redshifts < z < 0.2 so that af g (L s ,z s ) ~ af s (L s ,0), 
where z s ~ 0.17. The conversion of crg g from real space to redshift space is determined 
by 

a! g (L s ,z s ) = crf g (L s , z s )\Jk\J3(L s , z s )], (12) 



where 



K\P] = 1 + |s + \p\ (13) 
is Kaiser's factor \§T\, and 

P(L s ,z s )*^f^ (14) 
b[L s ,z s ) 

with 

+ z) 3 

n ^ = o n i Talo • ( 15 ) 
u m {l + zy + \ l A 

Furthermore, we assume that the mass fluctuations grow as a 8m (z) = a 8m (0)D(z), where 
D(z) is the linear growth rate (normalized to 1 at z — 0). As pointed out earlier, the 
linear growth rate is actually scale dependent in models with massive neutrinos. We 
have checked that this scale dependence is weak for the parameter range we consider, 
and evaluated D(z) at k — k, where k is the mean value of k for a spherical top-hat 
window function. With these assumptions we have 

b(L 8 ,0) 

b{L s ,z s ) = ——— . 16 
D{z s ) 

Given b(L s ,0), we can translate a given normalization crg m (0) of the matter power 
spectrum to cxf g (L s , z s ) for the galaxy power spectrum. 

We now carry out the following simple analysis: we fix n = 0.99, VL^ti 1 = 0.024 
(all values taken from table 1 in pQ for WMAP alone), and fit lo v , u m , h, a% m and 
b(L s ,0) to the combination of the 2dFGRS power spectrum data and the constraint 
c"8m^m 6 = 0.44 ± 0.10 from table 2 in |Tj. Furthermore, we add Gaussian priors 
Vl m h 2 = 0.14 ± 0.02, h = 0.72 ± 0.05 from WMAP, and look at the limit on 
m u,tot = 94:Q u h 2 eV for the cases of with and without a prior on b = b(L s ,0). The 
results reveal that the prior on uj m is crucial and the prior on h very important, but also 
that, at least in this simplified analysis, the effect of adding a prior on the bias is very 
small, as can be seen from figure ^2 Without a prior on uj m , no non-trivial limit on 
m v,tot is obtained. With just the cu m prior, the 95 % confidence limit is m u to i < 1.6 eV, 
and this improves to 1.1 eV when the prior on h is added. Adding the prior on b 
does not change the limit on m^ t ot, it stays at 1.1 eV. From the contour plots in 
figure H21 the degeneracy between uj v and b is small, especially when the prior on u m is 
included. This is not in contradiction to the analysis in |2Hj, because our treatment of 



Upper limits on neutrino masses from the 2dFGRS and WMAP: the role of priors 19 




Figure 11. Probability distributions for the total neutrino mass in our simplified 
analysis after marginalizing over u) m , h, and b for the cases of a prior on w m only, and 
priors on uj m and h. When a Gaussian prior b = 1.04 ± 0.11 is added, the resulting 
curve is indistinguishable from the dashed line. 
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Figure 12. Confidence contours (68 and 95 plane after marginalizing over w m and h 
for the cases of a prior on w m , and priors on both uj m and h. 



biasing is different: we have treated the redshift-space distortions explicitly, and then 
the constraint on o"g and uj m from the CMB breaks much of the degeneracy between 
uj v and b. This is because the redshift distortion itself depends on Q m , see equations 

(unm. 

In this simple analysis we get a 95 % confidence limit of m Uyto t < 1.1 eV. This is still 
some way from the WMAP limit of 0.71 eV, even with our very restricted parameter 
space, but consistent with the analysis in The WMAP analysis also used data from 
ACBAR and CBI [TH ITH] , and included the Lyman a forest power spectrum. The linear 
matter power spectrum inferred from the Lyman a fores probes smaller scales than the 
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2dFGRS and therefore has considerable power in constraining neutrino masses. We 
have seen that the most severe degeneracies of io v are with u m , n, h, and u^. The most 
serious one is with u m : without a prior on the physical matter density, one cannot get 
a non-trivial upper bound on m u ^ tot . This makes sense, as the matter power spectrum 
depends on f v = Q u /Q m and this is why the analysis was carried out in terms of this 
parameter in [TJ]. So the fact that WMAP (restricted to flat ACDM models and h > 0.5) 
provides a tight constraint u> m = 0.14 ±0.02 is crucial for constraining neutrino masses. 
WMAP also constrains the spectral index to a narrow interval around n = 1, and, just 
as importantly, by constraining h to be around 0.7 eliminates the possibility for MDM 
models with large neutrino masses to give good fits to the 2dFGRS power spectrum. 
The importance of the prior on h was also noted in [2*3] . 

4. Other cosmological probes of neutrino masses 

Direct probes of the total matter distribution avoid the issue of biasing and are therefore 
ideally suited for providing limits on the neutrino masses. Several ideas for how this can 
be done exist. In [HS] the normalization of the matter power spectrum on large scales 
derived from COBE was combined with constraints on as from cluster abundances and 
a constraint m u tot < 2.7 eV obtained, although with a fairly restricted parameter space. 
However, a 8 is probably one of the most debated numbers in cosmology at the moment 
|4*H] , and a better understanding of systematic uncertainties connected with the various 
methods for extracting it from observations is needed before this method can provide 
useful constraints. The potential of this method to push the value of the mass limit down 
also depends on the actual value of a 8 : the higher o~ 8 turns out to be, the less room 
there will be for massive neutrinos. The evolution of cluster abundance with redshift 
may provide further constraints on neutrino masses |4fij . The Lyman a forest provides 
constraints on the matter power spectrum on small scales, where the effect of massive 
neutrinos is most visible, and it was used in (HH] to derive a limit m Uttot < 5.5 eV, 
and it clearly played a role in the WMAP limit also. How to use this probe correctly 
in cosmological parameter estimation is, however, still being discussed jlHj. Massive 
neutrinos also suppress peculiar velocities on scales smaller than 50 h~ l Mpc, where 
they can be measured more accurately to nearby galaxies, however, non-linear effects 
on small scales causes complications. Finally, deep and wide weak lensing surveys will 
in the future make it possible to do weak lensing tomography of the matter density field 
[SHE], and in [Z2] it has been shown that one can probe neutrino masses below 0.1 eV 
in this way. However, this is under the assumption that the equation of state of the 
dark energy is known. 

5. Conclusions 

We have reviewed how the constraint on the neutrino mass in ^2] was obtained, and the 
recent improved limit from WMAP, paying attention to issues of priors and parameter 
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degeneracies. We have seen that one can derive fairly tight constraints on neutrino 
masses from the 2dFGRS power spectrum, provided that one has good constraints on 
u m , n, h, and u>b from independent data sets. 

We found that external constraints on the Hubble parameter, for example the HST 
Key project, is important in order to get a strong limit on the neutrino mass since, 
intriguingly, the out-of-fashion Mixed Dark Matter (MDM) model can give a reasonable 
description of the CMB and 2dFGRS data with Q u = 0.2, Q m = 1 and no cosmological 
constant if we allow low values of the Hubble parameter, h < 0.5. The importance 
of having a prior on h was also noted in j2H|. We note that the above MDM model 
is inconsistent with other cosmic measurements such as Supernovae Type la, baryon 
fraction in clusters, and the evolution of cluster abundance with redshift, so adding any 
one of these datasets to the analysis would eliminate MDM and improve the neutrino 
mass limit. 

We also considered the effect of the possible bias of the galaxy distribution with 
respect to the mass distribution on the neutrino mass limit. A scale-dependent bias 
has serious implications for the m Ujtot constraint, but based on semi- analytic galaxy 
formation models [H2| the scale-dependence of the bias is expected to be too weak to 
be of any major concern on the large scales used in the analysis of the 2dFGRS power 
spectrum. When the effects of redshift-space distortions on the 2dFGRS power spectrum 
are included in the analysis, there is almost no degeneracy between a constant, scale- 
independent bias factor and the neutrino mass. However, in our restricted analysis we 
did not get as good a neutrino mass constraint with 2dFGRS + WMAP priors as in the 
full analysis in [1. which suggests that the Lyman a forest power spectrum plays a role 
in pushing the constraint on m„ itot below 1 eV. 
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